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A search for a very-high-energy (VHE; > 100 GeV) 7-ray signal from self-annihilating particle 
Dark Matter (DM) is performed towards a region of projected distance r ~ 45 — 150 pc from 
the Galactic Center. The background-subtracted 7-ray spectrum measured with the High Energy 
Stereoscopic System (H.E.S.S.) 7-ray instrument in the energy range between 300 GeV and 30 TeV 
shows no hint of a residual 7-ray flux. Assuming conventional Navarro-Frenk- White (NFW) and 
Einasto density profiles, limits are derived on the velocity-weighted annihilation cross section (av) 
as a function of the DM particle mass. These are among the best reported so far for this energy 
range. In particular, for the DM particle mass of ~ 1 TeV, values for (ov) above 3 x 10 -25 cm 3 s _1 
are excluded for the Einasto density profile. The limits derived here differ much less for the chosen 
density profile parametrizations, as opposed to limits from 7-ray observations of dwarf galaxies or 
the very center of the Milky Way, where the discrepancy is significantly larger. 



INTRODUCTION 

The existence of particle DM is a widely accepted as- 
trophysical concept that is used to explain, e.g., the for- 
mation of large scale structure during the evolution of 
the universe. In the present universe, galaxies, among 
others, are believed to be embedded in DM halos, which 
have density profiles that are much more extended than 
the profiles of visible matter and also exhibit the highest 
DM density in their centers (e.g. [3-0|). This fact makes 
galactic centers promising targets to search for signals 
from annihilation or decay of hypothetical DM particles. 
In simulations of cosmological structure formation, cold 
DM, which, among others, might consist of non-baryonic, 
weakly interacting massive particles (m x ~ 10 GeV up 
to a few TeV) [4j , is able to reproduce the observed large 
scale structure of the universe [l|, 0] . They are expected 
to either annihilate or decay into standard model parti- 
cles, producing, among other particles, photons in the 
final state. These photons are predicted to exhibit a 



continuous spectrum and to cover a broad energy range 
up to the DM mass, with possibly narrow lines from di- 
rect XX ~~ ^ Jl/jZ annihilations or other spectral features 
overlaid (see, e.g., Therefore, 7-ray observations, 

in particular at very high energies where DM masses of 
a few 100 GeV and above are probed, are a promising 
way to detect DM in space and to constrain its particle 
physics properties. 

A successful DM observation strategy has to avoid 
sky regions with strong astrophysical 7-ray signals, and 
should focus at the same time on regions with an ex- 
pectedly large DM density. Besides dwarf galaxies, from 
which limits on (av) have been derived recently from ob- 
servations at high (100 MeV-20 GeV) [1, i| and very 
high energies da-El, the Galactic Center (GC) region is 
a prime target for DM search, both because of its prox- 
imity and its predicted large DM concentration (see, e.g., 
3]). As opposed to dwarf galaxies, however, the search 
for DM induced 7-rays in the GC is hampered by a strong 
astrophysical background. Especially at the very center, 
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there is the compact 7-ray source HESS J1745-290, co- 
incident with the position of the supermassive black hole 
Sgr A* and a nearby pulsar wind nebula jl6l4l8j . which 
makes the detection of 7-rays from DM annihilation diffi- 
cult Furthermore, there is a band of diffuse emission 
along the GC ridge [20] , which has already been used to 
derive upper limits on a DM induced signal 2ll423j. 

The 7-ray flux expected from annihilations of DM par- 
ticles of mass m x is given by the product of the integral J 
of the squared DM mass density p along the line-of-sight, 
the velocity- weighted annihilation cross section (av) , and 
the average 7-ray spectrum dN/dE^ produced in a sin- 
gle annihilation event [2j|. For a given DM mass and 
given branching ratios for annihilation into SM particles, 
the shape of the expected 7-ray spectrum, determined 
by dN/dE 1 , can be rather accurately computed, since 
only SM particle physics processes are involved. Measur- 
ing the velocity- weighted annihilation cross section (av), 
however, suffers from large systematic uncertainties on 
the astrophysical factor J. This is because in most indi- 
rect searches the main contribution of the total DM sig- 
nal arises from the very central region of the object under 
study, where the DM density profile peaks. In these inner 
regions, however, the profile is so far only poorly known. 
In particular, for a Milky- Way sized DM halo, the radial 
DM density profiles obtained by the Aquarius [l| and Via 
Lactea II [2| simulations can be described by Einasto and 
NFW parametrizations, respectively 25]. These profiles 
are shown in Fig. Q] as a function of galactocentric dis- 
tance r. Large differences between both parametrization 
s occur if they are extrapolated down to the very cen- 
ter of the halo, where the NFW profile is more strongly 
peaked. At distances > 10 pc, however, the difference is 
merely a factor of two, allowing one to put limits on (av) 
which do not depend strongly between either of the two 
parametrizations. 

Here we exploit this fact by searching for a VHE 7- 
ray signal from DM annihilation in our own Galaxy, 
in a region with a projected galactocentric distance of 
45 pc — 150_pc (corresponding to an angular distance of 
0.3° — 1.0° Jj, excluding the Galactic plane. In this way, 
contamination from 7-ray sources in the region is natu- 
rally avoided as well. 



METHODOLOGY 

The analysis is carried out using 112 h (live time) of 
GC observations with the H.E.S.S. VHE 7-ray instru- 
ment (see [271 and references therein) taken during the 
years 2004-2008. For minimum energy threshold, only 
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1 Here and in the following a distance of the GC to the observer 
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FIG. 1. Comparison of the Galactic DM halo profiles used 
in this analysis. The parameters for the NFW and Einasto 
profiles are taken from [25[. An isothermal profile |24|. ex- 
hibiting a flat DM density out to a galactocentric distance 
of 1 kpc, is shown for comparison. All profiles are normal- 
ized to the local DM density (p = 0.39 GeV/cm 3 [1| at a 
distance of 8.5 kpc from the GC). The source region and the 
region used for background estimation are indicated. Note 
that the predicted DM density is always larger in the source 
region, except for the isothermal profile, which is included for 
completeness. 



observations with zenith angles smaller than 30° are con- 
sidered. The mean zenith angle is 14°. To avoid possible 
systematic effects, pointing positions were chosen fairly 
symmetric with regard to the Galactic plane. The mean 
distance between the pointing position and the GC is 
0.7°, with a maximum of 1.5 c irc. Events passing H.E.S.S. 
standard cuts defined in 27J are selected for analysis. To 
minimize systematic uncertainties due to reduced 7-ray 
efficiency at the edges of the ~ 5° diameter field-of-view 
(FoV) of H.E.S.S., only events reconstructed within the 
central 4° are considered. The effective 7-ray collection 
area of the chosen event selection is « 1.7 X 10 5 m 2 at 1 
TeV. The total effective exposure of the utilized dataset 
at 1 TeV amounts to « 2.4 x 10 7 m 2 sr s. 

7-rays from DM annihilations are searched for in a cir- 
cular source region of radius R on = 1.0° centered at the 
GC. Contamination of the DM signal by local astrophys- 
ical 7-ray sources is excluded by restricting the analy- 
sis to Galactic latitudes |6| > 0.3°, effectively cutting 
the source region into two segments above and below the 
Galactic plane (see Fig. [2]). Simulations show that both 
for Einasto and NFW parametrizations the sensitivity 
varies only within a few percent when varying the source 
region size in the range 0.8° < R on < 1.2°. 

For ground-based VHE 7-ray instruments like 
H.E.S.S., any 7-ray signal is accompanied by a sizable 
number of cosmic-ray induced background events, which 
are subtracted from the source region using control re- 
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gions in the FoV of the observation. Fig. [2] visualizes 
details of the method, which is an evolution of the stan- 
dard reflected background technique [28[ adjusted for this 
particular analysis. By construction, background regions 
are located further away from the GC than the source 
region. This is an important aspect, since, unavoidably, 
a certain amount of DM annihilation events would be 
recorded in the background regions, too, reducing a po- 
tential excess signal obtained in the source region. For 
the NFW and Einasto profiles, the expected DM annihi- 
lation flux is thus smaller in the background regions than 
in the source region (cf. Fig.[lJ, making the measurement 
of a residual annihilation flux possible. Note, however, 
that for an isothermal halo profile, the signal would be 
completely subtracted. As far as the background from 
Galactic diffuse emission is concerned, its predicted flux 
[isij is significantly below the current analysis sensitivity, 
thus its contribution is not further considered in the anal- 
ysis. In any case, since its intensity is believed to drop 
as a function of Galactic latitude, 7-rays from Galactic 
diffuse emission would be part of a potential signal, and 
therefore lead to more conservative results for the upper 
limits derived in this analysis. 



RESULTS 

Using zenith angle-, energy- and offset-dependent ef- 
fective collection areas from 7-ray simulations, flux spec- 
tra shown in Fig. [3] are calculated from the number of 
events recorded in the source and background regions^. It 
should be stressed that these spectra consist of 7-ray- like 
cosmic-ray background events. Both source and back- 
ground spectra agree well within the errors, resulting in 
a null measurement for a potential DM annihilation sig- 
nal, from which upper limits on (av) can be determined. 

The mean astrophysical factors J src and Jb g are calcu- 
lated for the source and background regions, respectively. 
The density profiles are normalized to the local DM den- 
sity po = 0.39 GeV/cm 3 26]. Assuming an Einasto pro- 
file, J src = 3142 x p\ x d E and J~b g = 1535 xp^xdv, where 
Pb = 0.3 GeV/cm 3 is the conventional value for the local 
DM density and g?e = 8.5 kpc the distance of Earth to 
the GC. For a NFW profile, J SIC = 1604 x p| x d E and 
Jbg = 697 x p| x d,E are obtained. This means that for an 
assumed Einasto (NFW) profile, background subtraction 
reduces the excess DM annihilation flux in the source re- 
gion by 49 % (43 %), which is taken into account in the 
upper limit calculation. 

Under the assumption that DM particles annihi- 
late into quark-antiquark pairs and using a generic 
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FIG. 2. Illustration of the cosmic ray background subtrac- 
tion technique for a single telescope pointing position (de- 
picted by the star). Note that this position is only one of the 
several different pointing positions of the dataset. The DM 
source region is the green area inside the black contours, cen- 
tered on the GC (black triangle). Yellow regions are excluded 
from the analysis because of contamination by astrophysical 
sources. Corresponding areas for background estimation (red 
regions) are constructed by rotating individual pixels of size 
0.02° x 0.02° of the source region around the pointing posi- 
tion by 90°, 180°, and 270°. This choice guarantees similar 
7-ray detection efficiency in both the source and background 
regions. As an example, pixels labeled 1 and 2 serve as back- 
ground control regions for pixel 0. Pixel 3 is not considered 
for background estimation because it is located in an excluded 
region. Pixels in the source region, for which no background 
pixels can be constructed, are not considered in the analysis 
for this particular pointing position and are left blank. 



parametrization for a continuum spectrum of 7-rays cre- 
ated during the subsequent hadronization 3tJ 31 1, limits 
on (av) as a function of the DM particle mass are cal- 
culated for both density profiles (see Fig. 0]). These 
limits are among the most sensitive so far at very high 
energies, and in particular are the best for the Einasto 
density profile, for which at ~ 1 TeV values for (av) 
above 3 x 10 -25 cm 3 s _1 are excluded. As expected from 



2 The background spectrum is reseated by the ratio of the areas 
covered by source and background regions (cf. also [2810 , 



the astrophysical factors, the limits for the Einasto pro- 
file are better by a factor of two compared to those for 
the NFW profile. Still, the current limits are one order 
of magnitude above the region of the parameter space 
where supersymmetric models provide a viable DM can- 
didate (see Fig. [4j. Apart from the assumed density 
parametrizations and the shape of the 7-ray annihilation 
spectrum, the limits can shift by 30% due to both the 
uncertainty on the absolute flux measurement 27J and 



the uncertainty of 15% on the absolute energy scale. For 
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FIG. 3. Top panel: Reconstructed differential flux Fg rc /B g , 
weighted with E 2 ' 7 for better visibility, obtained for the source 
and background regions as defined in the text. The units are 
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Due to an energy-dependent selection 



efficiency and the use of effective areas obtained from 7-ray 
simulations, the reconstructed spectra are modified compared 
to the cosmic-ray power-law spectrum measured on Earth. 
Bottom panel: Flux residua F lea / AF lea , where F lea = fsrc — 
Fb s and AF ICS is the statistical error on F las . The residual 
flux is compatible with a null measurement. Comparable null 
residuals are obtained when varying the radius of the source 
region, subdividing the data set into different time periods 
or observation positions, or analyzing each half of the source 
region separately. 
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FIG. 4. Upper limits (at 95% CL) on the velocity-weighted 
annihilation cross-section (av) as a function of the DM par- 
ticle mass m x for the Einasto and NFW density profiles. 
The best sensitivity is achieved at m x ~ 1 TeV. For com- 
parison, the best limits derived from observations of dwarf 
galaxies at very high energies, i.e. Sgr Dwarf [10J, Will- 
man 1, Ursa Minor plj and Draco using in all cases 
NFW shaped DM profiles, are shown. Similar to source re- 
gion of the current analysis, dwarf galaxies are objects free 
of astrophysical background sources. The green points rep- 
resent DarkSUSY models [32[, which are in agreement with 
WMAP and collider constraints and were obtained with a 
random scan of the mSUGRA parameter space using the 
following parameter ranges: 10 GeV < Mo < 1000 GeV, 
10 GeV < Mi /2 < 1000 GeV, A = 0, < tan/3 < 60, 
sgn(M) = ±1- 



the latter case, apart from a displacement with regard to 
the DM particle mass scale, the limits shift up (down) if 
the 7-ray energy is overall under(over)estimated. 



SUMMARY 

A search for a VHE 7-ray signal from DM annihilations 
was conducted using H.E.S.S. data from the GC region. 
A circular region of radius 1° centered at the GC was cho- 
sen for the search, and contamination by astrophysical 
7-ray sources along the Galactic plane was excluded. An 
optimized background subtraction technique was devel- 
oped and applied to extract the 7-ray spectrum from the 
source region. The analysis resulted in the determination 
of stringent upper limits on the velocity- weighted DM an- 
nihilation cross-section (av) , being among the best so far 
at very high energies. At the same time, the limits do not 
differ strongly between NFW and Einasto parametriza- 
tions of the DM density profile of the Milky- Way. 
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